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ff^ ', 1. Introduction 

o 

^vq . One of the most important reeent discoveries in particlephysics is the ob serv ation 

of neutrino oscillations in atmospheric!^, solar^, reactor^ and accelerator l^l^J neu- 
trino experiments. Neutrino oscillations is a quantum-mechanical consequence of 
the neutrino mixing relation 



H 
.^.' J^il{x) =^UHVtL{x) {l = e,y.,T). (1) 

Here i'i{x) is the field of neutrinos with mass m-i, U is the 3x3 unitary PMN3 I I I 
mixing matrix. The left-handed flavor field I'lhix) enters into the standard leptonic 
charged current 

jf'^(a;)=2 Y. ^il{x)Jc.Il{x) (2) 

and determines the notion of a left-handed fiavor neutrino i>i which is produced in 
CC weak processes together with a lepton l^. The fiavor neutrino i^i is described 
by the mixed state 

3 

i=l 
1 
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where \ui) is the state of a neutrino with mass rrii and a definite momentum. 

The probabihty of the transition ui — >■ vi in vacuum is given by the standard 
expression (see Ref!^ 

2 
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iEit , 
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n'} = 



u, 



^ki 
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I'l 






i\u, 



(4) 



the neutrino source, and E is the neutrino energy. 

In the standard parameterization, the 3x3 PMNS mixing matrix is characterized 
by three mixing angles, ^12, ^23 and fJis, by a Dirac CP-violating phase 5 and by 
two possible Majorana CP-violating phases A2 and A3: 



U = I -S12C23 
S12S23 ■ 



C12C13 

- Ci2S23Sl3e*'^ 
- Ci2C23Si3e"^ 



C12C23 - 
-C12S23 



S12C13 

- Si2S23Sl3e 
Sl2C23Sl3e 



si3e 



-iS\ 



iS 



iS 



S23C13 
C23C13 



i^(A2,A3), (5) 



where 



Cab 



cos'dab and Sat = sin-dab- The diagonal matrix I?(A2,A3) 



diag(l, e^^, e^^) is present only if massive neutrinos are Majorana particles. The 
Majorana phases have an effect in processes which are allowed only if massive neu- 
trinos are Majorana particles and are characterized by a violation of the total lepton 
number, as neutrinoless double-beta decay (see Section HJ. Since neutrino oscilla- 
tions are flavor transitions without violation of the total lepton number, they do 
not depend on the Majorana phases . The neutrino oscillation probabil- 

ities depend only on the four mixing parameters i?i2, 'i?23 7 ^13 and S, and on two 
independent mass-squared differences Amj2 and Arn^^. From the analysis of the 
experimental data it follows that 



'12 



30 



lAni 



231- 



(6) 



In the case of three- neutrino mixing assumed in Eqs. ([T]) and ([3]), two neutrino mass 
spectra are possible: 



(1) Normal spectrum (NS) 



mi < 7712 < tn^] A771I2 ^ ^^23- 



(2) Inverted spectrum (IS) 



77i3 < rni < TO2; A?Ti]^2 "^ I Am 



131- 



(7) 



(8) 



The existing experimental data do not allow to establish, what type of neutrino 
mass spectrum is realized in nature. 

Let us introduce the "solar" and "atmospheric" mass-squared differences Am.^ 
and Am,^, respectively. For both spectra we have Atti^j = Am^. For normal (in- 
verted) spectrum we have Am,23 = Am^ (jAm^gl = Atti^). 
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From a three-neutrino analysis of the Super-Kamiokande data^, the values of 
the neutrino oscihation parameters in the case of a normal (inverted) mass spectrum 
are, at 90% C.L., 

1.9 (1.7) • 10-'"^ eV^ < Ami < 2-6 (2.7) • 10"^ eV^ 

0.407 < sin^ i?23 < 0.583, sin^ i^ig < 0.04 (0.09). (9) 

The results of the Super-Kamiokande atmospheric neutrino experiment have been 
fully confirmed by the long-baseline accelerator neutrino experiments K2K and 
MINOP. 

From the two-neutrino analysis of the MINOS i^^ — > f^ data, for the parameters 
Ami ^^d sin^ 2'i923 the following values were obtained: 

Ami = (2-32lg;o8) • 10"^ eV^ sin^ 2z?23 > 0.90. (10) 

From the combined three-neutrino analysis of all solar neutrino data and the 
data of the reactor KamLAND experiment, it was found that ■ J 

Ami == (7-50lj;-2o) ■ 10"^ eV^ tan^ ^^^^ = 0.452;°°i^^, sin^ ^^g = 0.020 ± 0.016. 

From a similar analysis performed by the SNO collaboration, it was obtained thalp^ 

Ami = (7.4ltl]?^) • 10-5 eV^ tan^ ^,, = 0.446;°°3°, sin^ ^13 = 0.025t°°i«. 

(12) 
The Daya Bay collaboration^^ measured recently with high precision the mixing 
angle ^13: 

sin^ ??i3 = 0.024 ± 0.004. (13) 

This is a 5.2ct evidence of a no n-ze ro value of ??i3 w hich confirms the previous 
measurements of T2KM^, MIN09^ and Double Chooa^. It also confirms earlier 
indications of a non-zero value of ^13 found in the analysis of the data of solar and 
other neutrino experiments (see Eqs. ([TT]) and ([12]) and Ref i— '-'^^ 20 21 2 2jj rpj^^ 
Daya B ay m easurement has important implications for theory^^ and experiment 
(see Ref.l^^. It opens promising perspectives for the observation of CP violation in 
the lepton sector and matter effects in long-baseline experiments, which could allow 
to determine the character of the neutrino mass spectrum. 

Several years ago an indication in favor of short-baseline P^ — > 9^ transitions was 
found in the LSND experiment '^. The LSND data can be explained by neutrino 
oscillations with 0.2 eV < Am^ < 2eV and 10"^ < sin^ 2^ < 4 • 10"^. Recently 
an additional {2a) indication in favor of short-baseline oscillations, compatible with 



the LSND result, was obtained in the MiniBooNE experiment!^. Moreover, the 
data obtained in old reactor short-baseline experiments can also be interpreted 
as indications infayor of oscillationa^ by using a new calculation of the reactor 
neutrino fluxe d ^i All these data (if confirmed) imply that the number of massive 
neutrinos is larger than three and in addition to the three flavor neutrinos J^e , ^'/j , i^r 
mixed sterile neutrinos Vsn ■■■ niust exist. 
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The problem of short-baseline neutrino oscillations and sterile neutrinos is a hot 
topic at the moment. Several new short-baseline reactor and a ccel erator experiments 
are aimed to check this possibility in the near future (see ref!^^. 

The absolute values of neutrino masses are currently unknown. The Maina^and 
TroitsH^ experiments on the high-precision measurement of the end-point part of 
the /3-spectrum of '^H decay found the 95% C.L. upper bounds 

mp < 2.3 eV (Mainz), m^ < 2.1 eV (Troitsk), (14) 

for the "average" neutrino mass (see Ref.^^^ 



From neutrino oscillation and tritium /3-decay data we conclude that 

(A) Neutrino masses are different from zero. 

(B) Neutrino masses are much smaller than the masses of charged leptons and 
quarks. 

(C) Neutrino masses are not (or not only) of Standard Model (SM) Higgs origin. 

Several mechanisms of neutrino mass generation have bee n proposed. It is widely 
believed that the most plausible one is the seesaw mechanisn i'-^ l ^^ l '^"!'-^ ' [ According 
to this mechanism, small neutrino masses are generated by new interactions beyond 
the SM which violates the total lepton number L at a scale much larger than the 
electroweak scale v = (V2Gf)"^/^ =i 246 GeV. 

If the seesaw mechanism is realized, the neutrinos Vi with definite masses are Ma- 
jorana particles and, consequently, the lepton number violating neutrinoless double- 
beta decay (Oi^/3/3-decay) of even-even nuclei, 

N{A,Z)^N{A,Z + 2) + e- +e-, (16) 

is allowed, where N(A^ Z) is a nucleus with nucleon number A and proton number Z. 
The knowledge of the nature of neutrinos with definite masses (Majorana or Dirac?) 
is extremely important for the understanding of the origin of small neutrino masses. 
Using large detector masses, high energy resolutions and low backgrounds, the ex- 
periments on the search for neutrinoless double-beta decay allow to reach unparal- 
leled sensitivities to extremely small effects due to t he Majorana neutrino m asses. 
In this brief review we consider this process (see also 



le Ma.iorana neutrino m ; 
1381391401411421431441451 ) 



2. Seesaw mechanism of neutrino mass generation 

In this Se ction we briefly discuss the standard seesaw mechanism of neutrino mass 
generatio n'^ I ^ . We consider a general approach based on the effective La- 

grangian formalisrrP^. Let us assume that the Standard Model is valid up to some 
scale A. If we include effects of physics beyond the SM, the total Lagrangian (in the 
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SM region) has the form 

£(A)=/:^^^ + ^-l04+„. (17) 

n>l 

The second term is a nonrenormahzable part of the Lagrangian. It is buih from SM 
fields and satisfies the requirement of SU{2) x f/(l) invariance. The operator O^+n 
has dimension M''"'""'. 

In the expansion (|17p of the non-renormahzable part of the Lagrangian in powers 
of 1/A, the most important term for neutrino physics is the first one, Cf^ — O5/A, 
which contains an operator of dimension five. This term can be buih from the 
leptons and Higgs doublets: 

^5'' = "X ^ I^^'lH] Yi'i [H^iLiLY] + h.c, (18) 

i',i 

for l^ /' = e, /i, T. Here 



l^lL 



H(+) 



are the lepton and Higgs doublets, H ~ iT2H* is the conjugated Higgs doublet, 
{LilY — C{Lil)^ is the (right-handed) charge-conjugated lepton doublet and Yin = 
Yui are dimensionless constants (presumably of order one). Here C is the charge- 
conjugation matrix (which satisfies the relations C7„ C~^ = — 7^ and C^ = —C). 

The Lagrangian (fT8)) does not conserve the total lepton number L. Let us stress 
that this is the only Lagrangian term with a dimension-five operator which can be 
built with the SM fields. 

The electroweak symmetry is spontaneously broken by the vacuum expectation 
value of the Higgs field 

From Eqs. (|18p and (j20p . we obtain the left-handed Majorana neutrino mass term 

C^^-lY.^i'LMi'^ii^iLT + h.c., (21) 

I'd 



where 



«2 



Mi% = -Yi,i (22) 

After the diagonalization of the symmetric matrix Y through the transformation 

Y = UyU^, U^U^l, y,k = y^S^k, (23) 

we obtain 

/:^--^^m,J?,t.,. (24) 
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Here 



and 






\y^, (25) 



uW+y^iuiyiLr. (26) 



I I 

From Eq. (|26p it follows that the field I'i satisfies the Majorana condition 

u, = ,.^^ CvJ. (27) 



Thus, Vi is the field of the Majorana neutrino with mass nii given by Eq. (j25p . 

From Eq. (pS]) . one can see that the flavor field i'il is connected to i^iL by the 
standard mixing relation 



J^lL 



= ^t/hi^^L, (28) 



where U is the unitary PMNS mixing matrix given in Eq. ([5]) in the standard 
parameterization, including the diagonal matrix of Majorana phases. 

The values of the neutrino masses are determined by the seesaw factor w^/A. 
Assuming that 7713 ~ 5 • 10^'^ eV (which is the largest neutrino mass in the case 
of a neutrino mass hierarchy mi <C m2 ^ ^13), we have A ~ 10^^ GeV. Thus, the 
standard seesaw mechanism of neutrino mass generation explains the smallness of 
neutrino masses by a violation of the total lepton number L in interactions due to 
physics beyond the SM at a very large (GUT) scale. 

The local effective Lagrangian (flS)) can be obtained by considering the possible 
existence of heavy Majorana leptons Ni with masses Mi ^ v, which are singlets of 
the SU{2)l X U{1)y gauge group of the SM. These heavy Majorana leptons can 
have the lepton number- violating Yukawa interaction with the standard lepton and 
Higgs doublets 

Cj = -V2Y,YuLilN,rH + h.c. (29) 

i,i 

At electroweak energies, the interaction (P^ generates the effective Lagrangian ([T5)) 
at second order of perturbation theory. We have 

E^-]^^^ = l>^- (30) 

i 

From this relation it follows that the masses Mj determine the scale of new physics. 

The seesaw mechanism based on the Lagrangian (|29l) is called "type I seesaw" . 
There are two other well-studied^ ways to generate the effective Lagrangian Dj^ 
and, consequently, the left-handed Majorana mass term (PT|) : through the interac- 
tion of the lepton and Higgs doublets with a heavy triplet scalar boson (type H 
seesaw) or with a heavy Majorana triplet fermion (type HI seesaw). 

Summarizing, if small neutrino masses are generated by the standard seesaw 
mechanism, we have the following consequences: 
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(1) Neutrinos with definite masses are truly neutral Majorana particles. 

(2) Neutrino masses are given by the seesaw relation (p5)) . Hence, the neutrino 
masses are suppressed with respect to the masses of charged leptons and quarks, 
which are proportional to v, by the small ratio v/ K. 

(3) The Majorana neutrino mass term is the only implication at the electroweak 
scale of a possible existence of heavy Majorana particles. 

(4) CP- violating decays of heavy Majorana particles in the early Universe could be 
the origin of the baryon asymmetry of the Universe (see Ref.!^^. 

3. On the theory of Oi^/J/J-decay 

In this Section we present a brief derivation of the matrix element of the neutrinoless 
double-beta decay process in Eq. (|16p . assumi ng that t his process is induced by 
Majorana neutrino masses and mixing (see Refs! ^ ^ *). 

The standard effective Hamiltonian of the process has the form 

Hi{x) = — =2eL(x)7ai^eL(a;)j"(x) +h.c. (31) 

Here Gp is the Fermi constant and j"(a;) is the hadronic charged current which 
does not change strangeness. In terms of the quark fields, the current j"(a;) has the 
form 

r{x) = 2cos^cUL{xh''dL{x). (32) 

The mixed flavor field Vehix) is given by the relation (P5|) with I = e: 

VeL (x) = V" Uei ViL {x) , (33) 



where U is the PMNS mixing matrix and Vi{x) is the field of the Majorana neutrino 
with mass m^, which satisfies the Majorana condition (|27l) . 

The process P^ is of second order in Gp, with the exchange of virtual neutrinos. 
The matrix element of the process is given by 

{f\S^\i) = -4 (^^ Np.Np, j d^x^d''x2Y,UL{pi)e'^''''lc.Ue^ 

x(0|r(j.,i(xi);.fjx2)|0)7jf/e.sJ(p2)e^''^^M^/|7^(^"(2;i)J^(a:2))|^.). (34) 

Here pi and p2 are electron four- momenta, J" {x) is the hadronic charged current 
in the Heisenberg representatiorO, Ni and Nf are the states of the initial and final 
nuclei with respective four- momenta Pi — {Ei,pi) and Pf = {Ef,pf), and NZ'^ = 
(27r)'^/^-\/2p° is the standard normalization factor. 

Tn Eq. I|34| l strong interactions are taken into account. 
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Fig. 1. Feynman diagram of the elementary particle transition which induces Ov/3/3-decay. 

Taking into account the Majorana condition (1271) . for the neutrino propagator 
we find the expressioio 



{0\T{ly,LiXl)P^Lix2))\0) 



(27r) 



d qe 



-iq{xi~X2) _ 



rrii 1-75 



^ ' q^ — m^ 2 



Performing the integration over x^, X2 and q'^ in Eqs. p4l) and psp . the matrix 
element of the process takes the form 

{f\S^\i) = 2z (^^ Np,Np,27r6{Ef + pI + P2 - E,)u{p,)j^^^{l + j5)Cu^{P2) 

X J d^x^d^X2e-^P^'^-^P-'^ J2 U!,m, J -03 - 



iq(xi—X2) 



J 

2-^ K_ -I- -nS -I- «" - K, -if 2^ 



En+pl + q^ - E,-ie 



En+p'i + q']-E,-ie 



(36) 



where q^ = y l*?!^ + ^ ^-nd E'n arc the energy levels of the intermediate nuclear 
state. 

This is an exact expression for the matrix element of Oi^/3/3-decay at second order 
of perturbation theory. In the following we consider major O^'' — > 0"*" transitions of 
even-even nuclei, for which the following standard approximations apply: 

(1) Effective Majorana mass approximation. 

Oi'/J/^-decay is due to the exchange of virtual neutrinos (see the diagram in 
Figll]). Taking into account that the average distance between nucleons in a 
nucleus is about 10"^'^ cm, the uncertainty relation implies that the average 



''The neutrino propagator is proportional to nii. This is connected to the fact that only left-handed 
neutrino fields enter into the Hamiltonian of weak interactions. Thus, in the case of massless 
neutrinos the matrix element of neutrinoless double /3-decay is equal to zero. This is a consequence 
of the g eneral theorem on the equivalence of the theories with massless Majorana and Dirac 
neutrino&^USlI 
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neutrino momentum is g ~ 100 MeV. On the other hand, from tritium exper- 
iments we have the upper bounds in Eq. (I14p . which constrain all the masses 
rrij to be smaller than about 2 eV. Therefore, the neutrino masses can be safely 



neglected in the denominators in Eq. (pS)) and we have g^ = \ \<t\^ + "^? — 9, 
with q=\<f\. 

Thus, from Eq. ([55| it follows that in the matrix element of {)v fi ji-Aec&j the 
neutrino properties and the nuclear properties are factorized and the neutrino 
masses and mixing enter into the matrix element in the form of the effective 
Majorana mass 

mpi3 = '^U^im,. (37) 

i 

(2) Long-wavelength approximation. 

We have \pkXk\ < \pk\R (k = 1,2), where R ~ 1.2^^/'^ • 10^^^^ cm is the radius 
of a nucleus with nucleon number A. Taking into account that |pfc | ^ 1 MeV, we 
have IpfeXfcl <C 1. Thus, we have e~*Pi^i~*P2^2 ^ i (this approximation means 
that electrons are produced in S'-states). 

(3) Closure approximation. 

The energy of the virtual neutrino, q ~ 100 MeV, is much larger than the 
excitation energy £"„ — Ei. Thus, the energy of the intermediate states En can 
be approximated by an average energy E. In this approximation, called "closure 
approximation" , we have 

{Nf\J'^{x,)\N„){Nn\JHx2))m ^ {Nj\J'^{x^)JP{x2))\N,) 
En+pl + q^-E,-ie E + pl + q-Ei^ie ' 

Taking into account these approximations and considering commuting hadronic 
currents (see Eqs. (|4ip and (|1^ below), for the matrix element of Oj^/3/3-decay we 
obtain the expression 

(/|5(2)|z) = Sttz f^) mpf,Np,Np,u{Pi){^ + l5)Cu^{P2) x 

d''xid^X2{Nf\J''{xi)Ki\xi - X2\)Mx2))\N,)6iEf+p1+p''^ - E,). 



Here 



piq{Si-x2) 
K{\xi-X2\) = —^ I d^q- 



(2^)^^ q E + q-iM, + Mf)/2 



(39) 



(40) 



where Mi{Mf) is the mass of the initial (final) nucleus. 

In the calculation of the hadronic part of the matrix element of Ovfi/S-decay, 
the following approximate expression for the effective charged current J°'{x) = 
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(jO(x),J(x))isuse(iSl 

A 

and 

A 



gA{q )Sn + 9M{q )i 



2N-0'n X 9 



2mr. 



- 9p{q^ 



, (ct„ • q) q 



2mr, 



(41) 



(42) 



Here, aJ^ and r^ are Pauli matrices acting, respectively, on the spin and isospin 
doublets of the n nucleon, tj^ — (ti + iT2)/2, r„ is the coordinate of the n nucleon, 
nip is the proton mass, gv{q^),gA{q^),gM{q^) and gp{q^) are the vector, axial, 
magnetic and pseudoscalar weak form factors of the nucleon. From the conserved 
vector current (CVC) and partially conserved axial current (PCAC) hypotheses, it 
follows that 



gv{q^)^Fl{q^)^F-{q^), gM{q') = Fi{q'yF^{q^), gp{q^) = 



2mpgA 

q^ + m^ 



, (43) 



where Ff^^' and Fj " are the Dirac and Pauli electromagnetic form factors of the 
proton (neutron) and g^i ~ 1.27 is the axial coupling constant of the nucleon. 

The expressions (|4ip and (|42l) can be obtained from the one-nucleon matrix 
element of the hadronic charged current. For the number density of nucleons in a 
nucleus, the following approximate expression is used: 



*(f )7°*(f) = ^ (5(a 



0- 



(44) 



n=l 



The nuclear matrix element (NME) M^^ , which is the integrated product of two 
hadronic charged currents and a neutrino propagator, is a sum of a Fermi (F), a 
Gamow- Teller (GT) and a tensor (T) term: 



k,l 



Hpirki) 



9\ 



+ HcTirki)^k ■ <?/ - HT{rki)S, 



kl 



\o'. 



(45) 



Here Ski — 3((Tfc • rki ) (<?i • ^ki ) ^ f?fc • (?; , with rki = rk — ri, and the neutrino potentials 
HF,GT,T{j'ki) are given by the expressions 



HF,GT,T{rkl) = -R / 
TT Jo 



jn,o,2iQ'^ki)hF,GT,T{q'^)q 
q + E-{M, + Mf)/2 



dq, 



(46) 



where R is the radius of the nucleus, and the functions hp^QT^riq^) are combinations 
of different form factorfjj. 



"^The functions hp_cT,T{<i ) can be found in Ref! 



El 
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Taking into account the Coulomb interaction of the electrons and the final nu- 
cleus, for the total width of Ot^/3/3-decay we find the general expression 

r""- = i - G'^iQ, Z) lAfOn' ^^, (47) 

where G^'^{Q, Z) is a known integral over the phase space, Q = Mi — Mf — 2 Trig is 
the Q-value of the process, and rUe is the electron mass. The numerical values of 
G^'''{Q, Z), Q and the natural abundance of several nuclei of experimental interest 
are presented in Table [T] 

4. Effective Majorana mass 

The effective Majorana mass m.^/3 is determined by the neutrino masses, the mixing 
angles and the Majorana phases. In this Section we discuss which are the possible 
values of the effective Majorana mass which can be obtained taking into account the 
information on the neutrino mass-squared differences and mixing angles obtained 
from neutrino oscillation data. 

In the standard parameterization ([5]) of the mixing matrix, we have 

|m^^| = jcos^ i9i2 cos^ 'dismi + e^'"^^ sin^ t9i2 cos^ ^131712 + e^'"^^ sin^ 1^13^3 1 , 

(48) 
where 0112 and ai3 are, respectively, the phase differences of Ue2 and Ues with 
respect to Uei'. ai2 = A2 and ai3 — X3 — 6 in the standard parameterization ([5]) of 
the mixing matrix. Therefore, Oi^/3/3-decay depends not only on the mixing angles 
and Dirac CP-violating phase, but also on the Majorana CP-violating phases. This 
is in agreement with the discussion after Eq. ([5]), since the total lepton number is 
violated in Oi^/3/3-decay. 

In the case of a NS, the neutrino masses 7712 and 7713 are connected with the 
lightest mass toi by the relations 



'TI2 = \/ "^1 + Arn^, m^ — y mf + Arn^ + Arn^. (49) 

On the other hand, in a IS m^ is the lightest mass and we have 



mi — y m,^ + Am"^, m,2 = y m,\ + Am^ + Ami. (50) 

Table 1. The values of G^'^ {Q, Z), Q and natural abunda nce o f the initial isotope for several /3/3-decay 
processes of experimental interest. Table adapted from Ref.l^^. 



/3/3-decay 


G"" 


Q 


nat. abund 


experiments 




[10-14 y-l] 


[keV] 


l%] 




«Ca ^ 48 Ti 


6.3 


4273.7 


0.187 


CANDLES 


76 Qg ^ 76 gg 


0.63 


2039.1 


7.8 


GERDA, Majorana 


82gg ^ 82Kr 


2.7 


2995.5 


9.2 


SuperNEMO, Lucifer 


IOOMo ^ iooru 


4.4 


3035.0 


9.6 


MOON, AMoRe 


ii6cd^ii6Sn 


4.6 


2809 


7.6 


Cobra 


iSOrpg ^ I30xe 


4.1 


2530.3 


34.5 


CUORE 


I36xe ^ l36Ba 


4.3 


2461.9 


8.9 


EXO, KamLAND-Zen, NEXT, XMASS 


iSONd ^ I50sm 


19.2 


3367.3 


5.6 


SNO-I-, DCBA/MTD 
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Fig. 2. Value of the effective Majorana mass |m^^| as a function of the hghtest neutrino mass 
in the normal (NS, with mrr,;„ = rm ) and inverted (IS, with mmin = "is) neutrino mass spectra 
before and after the Daya Baj measurement of i9i3 in Eq. II13I I. The curr ent upper bound on 
l^i/s^l (see Eqs. I|70ll . II72II and lj74p 'l and the cosmological bound (see Ref!^^ on ^^rm ~ Smmin 
in the quasi-degenerate region are indicated. 



Figure [2] shows the val ue of the effective Majorana mass |m^^| as a function of 
the hghtest neutrino mas a I I in the normal and inverted neutrino mass spectra 
before and after the Daya Bay^ measurement oi'dis in Eq. P^ . We used the values 
of the neutrino oscillation parameters obtained in the global analysis presented in 
ReM 



Atoi2 



_ f.q + (0.20,0.40,0.60) ,„- 
'•°^-(0.18,0.35,0.50) ^ ^" 



'eV^ sin^z9i2 =0.312 



+ (0.017,0.038,0.058) 
-(0.015,0.032,0.042)' 



and in the NS 

ATOj3 — ^-^"-(0.16,0.25,0.36) 

whereas in the IS 



2-9 c;n+(009.0-18'0-26) , „- 



•eY' 



sin^z? 



13 



0.013 



+ (0.007,0.015,0.022) 
-(0.005,0.009,0.012)' 



^13 



-(0.08,0.18.0.27) -,n-3 Tr2 ■ 2 n nnicH 

-(0.09,0.17:0.27) ^10 cV , Sm ?9i3 = 0.016_ 



(0.008,0.015,0.023) 
(0.006,0.011,0.015)- 



(51) 



(52) 



(53) 



The three levels of uncertainties correspond to {la, 2a, 3a). In the "After Daya Bay" 
plot in Fig [5] we replaced the value of i?i3 in Eqs. ([5^ and ([55)) with that measured 
by the Daya Bay Collaboration in Eq. (fT5)) . The uncertainties for \mi3p\ have been 
calculated using the standard method of propagation of uncorrelated errors, taking 
into account the asymmetric uncertainties in Eqs. (|5ip - (l53l) . 

In the following we discuss the predictions for the effective Majorana mass in 
three cases with characteristic neutrino mass spectra: 



(1) Hierarchy of neutrino masseo 



mi ^ m2 ^ TO3. 



(54) 



Quarks and charged leptons have this type of mass spectrum. 
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(2) Inverted hierarchy of neutrino masses: 

ma < mi < m2. (55) 

(3) Quasi-degenerate neutrino mass spectrum: 

/a — 2^ imi<m2<m3 (NS), 

\/Ami < mo ~ <^ , , 56 

[m3<mi<m2 (IS), 

where m,o is the absohite mass scale common to the three masses. As one can see 
from Fig. [2l the Daya Bay measurement of dis has a visible impact on the value of 
|m^^| only in the case of a hierarchy of neutrino masses, discussed in the following, 
because only in that case the contribution of the largest mass m^ , which is weighted 
by sin^ i?i3, is decisive. 

4.1. Hierarchy of neutrino masses 

In this case we have 

mi ^ vAmf , m2 ~ y/Am^, m^ ~ \J Am\. (57) 

Thus, m2 and 7713 are determined by the solar and atmospheric neutrino mass- 
squared differences. Neglecting the contribution of mi to the effective Majorana 
mass, from Eq. (|48|) we find 



l™/3/3| 



sin^ i?i2 v/Am2 ^ e^"^'-' sin^ di-i VArn^ , (58) 



where 023 is the phase difference between Ues and Ue2'- 0^23 = 0^13 — cn2 = A3— (5 — A2 
in the standard parameterization ([5]) of the mixing matrix. 

The first term in Eq. (|58| is small because of the smallness of Am^ . On the other 
hand, the contribution of the "large" Attl^ is suppressed by the small factor sin^ -di^. 
Hence, both terms must be taken into account and cancellations are possible, as 
shown in Fig. [H 

As one can see from Fig. [21 in the case of a hierarchy of neutrino masses we have 
the upper bound 

lm/3/3 1 < sin^ 79i2 v/Am2 + sin^ 1^13 \/Am2 < 5 • 10"^ eV, (59) 

which is significantly smaller than the expected sensitivity of the future experiments 
on the search for Oi//3/3-decay (see Section[6|). This bound corresponds to the case of 
g2ia23 = 1. It is slightly increased by the Daya Bay measurement of -di^ in Eq. I\l'6\i . 
because the additive contribution of sin^ 1^13 -^/Am^ in Eq. (f58|) is increased. On 
the other hand, one can see from Fig. [2] that the lower bound on \nipp\ for mi <C 
10"'^ eV, which corresponds to e^*"^^ = —1, is slightly decreased by the Daya Bay 
measurement of i9i3, because the increased contribution of sin^ ■^l3^JArr>^ in this 
case is subtracted. 

From Fig. [2] one can also see that when the contribution of m-i is not negligible, 
there can be cancellations among the three mass contributions. The two extreme 
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cases in which cancellations can happen are the following ones in which CP is 
conserved: 

g2iQi2 _ _2 and e^*"^^ = +1. The value of mi for which cancellations suppress 
|to/3^| is slightly decreased by the Daya Bay measurement of z9i3, because the 
larger value of sin fJisTOa adds to the contribution of nii . Hence, a smaller value 
of mi is required to cancel the sum of the contributions of mi and 7713 with the 
opposite contribution of m2 . 

g2iQi2 _ _2 aj^fj g2jQi3 _ _2. The value of mi for which cancellations suppress 
Imppl is slightly increased by the Daya Bay measurement of 1^13, because the 
larger value of sin^ ^137713 adds to the contribution of m2- Hence, a larger value 
of TOi is required to cancel the contribution of mi with the opposite sum of 
contributions of m,2 and 7713. 

Figure [2] shows thatO the two effects lead to a slight widening of the cancellation 
band after the Daya Bay measurement of i9i3. 

4.2. Inverted hierarchy of the neutrino masses 

In this case, for the neutrino masses we have 

7713 < vArof , 7771 — VA7772, 7772 ~ V^^ 1 + r. A % — \/ ^"^\- (60) 

V 2A7n^ ) 

In the expression of |777^^|, the contribution of the term 7773 sin 7^13 can be safely 
neglected. Neglecting also the small contribution of sin^ 7^13, from Eq. (|48|) we find 



|777^^ I ~ ^ Am2 (1 - sin^ 27?i2 sin^ aia) . (61) 

The phase ai2 is the only unknown parameter in the expression for the effective 
Major ana mass in the case of a inverted mass hierarchy. 
From Eq. (|6ip we find the following range for |777^^|: 

cos279i2\/Ato2 < \msifi\ < ^J l\ml. (62) 

The upper and lower bounds of this inequality correspond to the case of 
CP-invariancc in the lepton sector. In fact, CP invariance implies that (see 
Refs.5(I3159) 

e2--=772 77t, (63) 

where 77^ — ±7 is the CP parity of the Majorana neutrino v^. If 772 = 771, we 
have ai2 = 0,7r (the upper bound in the inequality (j62p ). If 772 = —771 we have 
ai2 — ±tt/2 (the lower bound in the inequality (p^ ). 

" We are very grateful to Michele Frigerio for pointing out a mistake in the cancellation band 
presented in the first arXiv version of this paper and in its published version (Mod. Phys. Lett. A 
27 (2012) 1230015). 
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From the existing neutrino oscillation data, we find the following range for the 
possible value of the effective Majorana mass: 

10"^ < |m/3^| <5-10"2cV. (64) 

The anticipated sensitivities to Ifn^^l of the future experiments on the search for 
the Oi^/3/3-decay are in the range (|64|) (see Section |6]). Thus, the future Oi^/3/3-decay 
experiments will probe the Majorana nature of neutrinos if a inverted hierarchy of 
neutrino masses is realized in nature. 

4.3. Quasi- degenerate neutrino mass spectrum 

Neglecting the small contribution of sin^ z9i3 in Eq. P5|) . in the case of a quasi- 
degenerate neutrino mass spectrum we obtain 



|m,3/3| ~ moy 1 - sin^ 2^12 sin^ai2, (65) 

where ttiq is the unknown absolute mass scale of neutrino masses (see Eq. (|56p) 
and ai2 is the phase difference between Ue2 and Uei- ai2 ~ ^2 in the standard 
parameterization ([5]) of the mixing matrix. Thus, in this case |?7i^^| depends on two 
unknown parameters: itiq and ai2. 

From Eq. ([65)) . we obtain the following range for the effective Majorana mass: 

cos2z?i2r7io < \mi3p\ < ttiq. (66) 

If Oz^/3/3-decay will be observed and the effective Majorana mass will turn out to 
be relatively large (|n^^/3| ^ i/ArriJ), we will have an evidence that neutrinos are 
Majorana particles and their mass spectrum is quasi-degenerate. In this case, we 
have 

|m^^|<mo< J^^^2.8|m^^|. (67) 

Information about the value of the mass scale will be inferred from the data of 
the future tritium /3-deca.y experiment KATRIN "" "^ and from future cosmologi- 
cal observations. The sensitivity of the KATRIN experiment to the neutrino mass 
scale is expected to be about 0.2 eV, which the same as the sensitivity to r/i^ in 
Eq. p5t . since in the quasi-degenerate case nijs ~ mo- Cosmological observations 
give information on the value of the sum of the neutrino masses ^^ ?77,i ~ Smg in the 
quasi-degenerate case. The existing cosmological data imply the bound ^^ rrii < 0.5 
eV (see Ref.l^^. It is expected that future cosmological observations will be sensitive 
to X^i '^j in the range (6 x 10"'^ — 10~^) eV (see, for example, Ref. ). 

5. Nuclear matrix elements 

The effective Majorana mass |?7i^^| is not a directly measurable quantity. The mea- 
surement of the half-life of Oi^/3/3-decay gives the product of the effective Majorana 
mass and the nuclear matrix element (see Eq. ()47p ). Hence, in order to determine 
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Fig. 3. Values of the NME calculated with the methods in Tab. [2] 



the effective Majorana mass one must calculate the nuclear matrix elements (NMEs) 
of Oj^/3/3-decay, which is a complicated nuclear many-body problem. Five different 
methods are used at present. In this short review we do not describe these methods 
and we do not discuss the advantages and disadvantages of each of them. We only 
present the references to the original papers in Tab.[2]and the latest results in Fig.[3l 
From Fig. [3] we reach the following conclusions: 

(1) The LSSM value of each NME is typically smaller than the corresponding one 
calculated with other approaches. Moreover, the LSSM value of each NME 
depends weakly on the nucleus, except for the double-magic nucleus ^^Ca. If 
Oi^/3/3-decay of different nuclei will be observed in future experiments, this char- 
acteristic feature of the LSSM can be checked, because the LSSM predicts the 
following ratio of half-lives of different nuclei: 



ri%(^i,Ai) 






A. 






(68) 



Table 2. Methods of calculation of nuclear matrix elements of Oi^/3/3-decay. 
Method References 



Quasi-particle Random Phase Approximation (QRPA) 

Energy Density Functional method (EDF) 

Projected Hartree-Fock-Bogoliubov approach (PHFB) 

Interacting Boson Model-2 (IBM-2) 

Large-Scale Shell Model (LSSM) 



.67.68. 

I69I70I 

I71I72I73I 
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(2) There is a large discrepancy between the values of NMEs calculated with dif- 
ferent approaches. The ratios of the maximal and minimal values of each NME 
are 3.1 (^^Ca), 2.4 (^^Ge), 2.0 (^^Se), 3.7 (^^Zr), 1.8 (IOOMo), 1.3 ("^Cd), 1.8 
(i24Sn), 1.9 (i28Te), 2.1 (""Te), 1.9 {^^^Xe), 2.3 (i^ONd). Therefore, the situa- 
tion with the calculation of the Oi^/3/3-decay NMEs is obviously not satisfactory 
at present. Further efforts and progress are definitely needed. 

6. Neutrinoless double-beta decay experiments 

Many experiments searched for neutrinoless double-beta decay without finding an 
uncontroversial positive evidence. The most stringent lower bounds on the half-lives 
of the decays of ^^Ge, -'^'^"Te and ^°°Mo have bee n obta ined, correspondingly, in the 
Heidelberg-MoscowJ^, CuoricincPS! and NEMOd^^^lSQI experiments. 

In the Heidelberg-Moscow experimentM-^ Germanium crystals with a 86% en- 
richment in the /3/3-decaying isotope ^^Ge were used. The total mass of ''^Ge was 
11 kg, with a low background of 0.11 counts/(kg keV y). After 13 years of running 
(with a 35.5 kg y exposure) no /3/3-peak at Q—2039 keV was found. The resulting 
half-live is 

T^J^C^Ge) > 1.9 X lO^^y (90%CL), (69) 

which implies tha1o 

Ito^^I < (0.22-0.64)eV. (70) 

In the cryogenic experiment Cuoricino-^ Te02 bolometers were used, with a 
total mass of 11.34 kg of ^'^°Te. The background was 0.17 counts/(kg keV y). After 
a 19.75 kg y exposure the following lower bound was obtained: 

T{';2("°Te) > 2.8 X lO^V (90%CL), (71) 

which corresponds to 

l"J/30l ^ (0.30-0.71)eV. (72) 

In the NEM03 experimenlM-^'^yi ^j^g cylindrical source was divided in sectors 
with enriched i°°Mo (6914 g), ^^Se (932 g) and other ;3/3-decaying isotopes. The 
two emitted electrons were detected in drift cells and plastic scintillator. No Oz^/3/3- 
decay was observed. The half-life of Oi^/3/3-decay of ^"°Mo have been bounded by 

T°;2('°"Mo)> 1.1 X 10^4 y (90%CL). (73) 

The corresponding limit for the effective Majorana mass is 

|w/3/3| < (0.44-1.00)eV. (74) 



Some participants of the Heidelberg-Moscow experiment claimedsiJthe observation of Oj^/3/3-decay 
of '"^Ge with half-life T°''C'^Ge) = (2.23to'3i) x 10^5 y (with 51.39 kg y exposure). From this 
result the authors found \mpp\ = 0.32 it 0.03 eV. This claim will be checked by the GERDA 
experiment" using the same Oj^/3/3-decaying nucleus. 
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Several new experiments on the search for Oi^/3/3-decay of different nuclei are 
currently running or in preparation. In the following we discuss briefly some of 
them (for more detailed presentations of future experiments see Ref ! I i) . 

In the GERDA experimenl!^, started in 2011, 18 kg of enriched germanium 
crystals (with 86% of the /3/3-decaying isotope ''^Ge) are used. The expected back- 
ground in the Phase-I of the experiment is 10"^ counts/(kg keV y). After one year 
of running it is expected to reach a sensitivity of T^Y^C^^Ge) — 2.5 x 10^^ y, which 

should allow to check the claim made in ReLl^. 

During the Phase-II of the GERDA experiment (expected to start in 2013), 
an array of enriched Germanium crystals (with 40 kg of ^^Ge) will be cooled 
and shielded by liquid Argon of very high radiopurity. A low background (10"'^ 
counts/ (kg keV y)) is expected. After 5 years of data taking, in the Phase-II of the 
experiment a sensitivity of T^Y^(^^Ge) ~ 1.9 • 10^^ y is expected. The corresponding 
sensitivity to the effective Majorana mass is Impnl ~ (7.3 • 10^^ — 2.0 • 10^^) cV. 

In the cryogenic CUORE experiment Te02 bolometers are used both as source 
and as detector. In the Phase-I of the experiment (started in the end of 2011) the 
target mass is 10.8 kg of ^■^°Te. In the Phase-II (expected to start in 2014) the 
target mass will be 206 kg of ^"^"Te. The expected background in this phase will be 
10~^ counts/(kg keV y). After 5 years of data taking a sensitivity of T^Y^{^^°Te) = 
1.6 • 10^^ y will be reached, which corresponds to |m^^| ~ (4.0 — 9.4) • 10~^ eV. 

In the KamLAND-Zen experimentl^^, the 0^/3/3-decay of ^'^^Xe will be studied. In 
this experiment enriched Xe (with 91% of the /3/3-decaying isotope -'^^^Xe) dissolved 
in liquid scintillator will be placed in a balloon (3.4 m in diameter) at the center 
of the KamLAND detector. In the first phase of the experiment (started in 2011), 
the source mass is 364 kg of ^"^^Xe. In the second phase (scheduled for 2013) 910 kg 
of -'^^^Xe will be utilized. After 5 years of data taking it will be possible to reach a 
sensitivity to |to/3^| in the region of the inverted hierarchy (|m^^| ~ 2.5 • 10~^ eV). 

In the running EXO experimenll^ the decay ^■^^Xe -> ^■^^Ba+e^+e^ is searched 
for. In the first phase of the experiment (EXO-200) the mass of the fiducial volume 
is about 150 kg of hquid Xenon enriched to 80.6% in the /3/3-decaying isotope ^'^^Xe. 
After two years of data taking a sensitivity |m^/3| ~ (8.7 • 10^^ — 2.2 • 10"-^) eV is 
planned to be achieved. The full EXO experiment will consist of about 1 ton of 
enriched liquid Xenon. With Ba^ tagging, a very low background of about 10~^ 
counts/ (kg keV y) will be reached. After 5 years of data taking, it is expected to 
reach a sensitivity of T^Y^l^^^Xe) ~ 10^^ y, which corresponds to |7tt.^/3| ~ (1.6 — 
4.0) • 10-2 eV. 



7. Conclusions 

If massive neutrinos are Majorana particles, neutrinoless double-beta decay of ™Ge, 
-'^™Mo, ^'^"Te, ^'^^Xe and other even-even nuclei is allowed. However, the expected 
probability of Oi^/3/3-decay is extremely small, because: 
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(1) It is a process of second order in the Fermi constant Gp- 

(2) Since the Hamihonian of weak interactions conserves hehcity, the amphtude of 
Oi^/3/3-decay is proportional to the very smaU factor 

^, (75) 

which comes from the neutrino propagator. Here m^^ = ^^ U'^^nii is the effec- 
tive Majorana mass (< 1 eV) and g is the average neutrino momentum (^ 100 

MeV). 

The expected half-hves of Oi^/?/3-decays depend on the decaying nucleus and are 
typically larger than 10^** — 10^^ years. Therefore, the observation of this rare process 
is a real challenge. 

The effective Majorana mass (and consequently the matrix element of the pro- 
cess) depends on the character of the neutrino mass spectrum. 

In the case of a quasi-degenerate spectrum, the expected value of to^^ is rela- 
tively large. This case is partly excluded by the data of the performed Oz^/3/3-decay 
experiments and by cosmological data (see Fig. [2]). It will be further explored by 
GERDA, KamLAND-Zen, EXO, CUORE and other experiments. 

In order to reach the region of the inverted neutrino mass hierarchy, with 10^^ ^ 
\mpij\ < 5- 10~^ eV, the construction of large detectors (^ 1 ton) and about 5 years 
of data taking will be required. 

We considered here the Oz^/3/3-decay induced by the standard mechanism of ex- 
change of light Majorana neutrinos between n-p-e~ vertices. From neutrino oscil- 
lation data it follows that if neutrino with definite masses are Majorana particles 
this decay mechanism is realized if there is no cancellation of the different mass 
contributions (as shown in Fig. [2l cancellations can happen in the normal scheme) . 

As discussed in Section [2l the neutrin o mass mech anism of Oi^/3/3-decay is pre- 
dicted by the standard seesaw mechanisn J'-^ l *--*'^ ! '^" ! '-^ ' [ However, additional sources 
of violation of the total lepton number L are possible (see Ref.l^ and references 
therein). If L is violated at the TeV scale these additional mechanisms could give 
contributions to the matrix elements of the 0://3/3-decay comparable with the con- 
tribution of the light Majorana neutrino mass mechanism. 

Let us consider as an example the violation of L due to R-parity violating in- 
teractions of SM and SUSY particles. In this case, Oz//3/3-decay is induced by the 
exchange of a heavy Majorana SUSY neutralino. The product of n-p-e~ vertices is 
given by the factor 

where A characterizes the scale of the masses of SUSY particles and mw is the mass 
of the W-boson. The factor (I76p must be compared with the corresponding factor 
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which appears in the case of the Majorana neutrino mass mechanism. Taking into 
account that q ~ 100 MeV and assuming that |m^^| ~ 10"^ eV, we come to the 
conchision that Eqs. (|76)) and ((77)) are comparable if A is of the order of a few TeV. 

If the Oi//3/3-decay of different nuclei will be observed in future experiments, it 
will be possible to probe the presence of different mechanisms which can generate 
the process. 

Finally, let us emphasize that the search for Oi//3/3-decay is a powerful practical 
way to solve one of the most fundamental problem of modern neutrino physics: are 
neutrinos with definite masses vi truly neutral Majorana particles or are they Dirac 
particles possessing a conserved total lepton number? 
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Note Added 

After the completion of this review, the EXO collaboration published in 
arXiv: 1205.5608 the important first result of EXO-200. With an exposure of 
32.5 kg y, they obtained T^J'^i^^^Xe) > 1.6 x lO^^y at 90% CL, corresponding 
to Impisl < (0.14- 0.38) eV. 
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